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Abstract 

The ordering transitions of a 2D lattice liquid characterized by a single favoured local structure 
(FLS) are studied using Monte Carlo simulations. All eight distinct geometries for the FLS are 
considered and we find a variety of ordering transitions - first order, continuous and multi-step 
transitions. Using a microcanonical representation of the freezing transition we resolve the dual 
influence of the local structure on the ordering transition, i.e. via the energy of the crystal and 
the entropy cost of structure in the liquid. 
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I. INTRODUCTION 



It is reasonable to presume that among the local structural arrangements that are found 
to be particularly stable in the liquids will be the ones found in the crystals. It is difficult to 
see how a liquid could freeze if this was not so. Let us imagine that the crystal is assembled 
so as to maximise the density of just a single local structure. It would follow that a minimal 
condition on any Hamiltonian used to model both the liquid and crystal phases is that it 
must stabilize this particular favoured local structure. The simplest such model would be 
the one that only stabilized this one local structure. This is a simpler picture than is usually 
invoked for liquid structure. There is, for example, an extended tradition 1 2 of focussing 
on non-crystalline favoured local structures in order to explain the kinetic stability of the 
supercooled liquid with respect to freezing. We note, in passing, that, while reasonable, there 
is no compelling evidence to suggest that the presence of a competing structure is necessary 
to explain liquid met ast ability. As we shall see, our simplest case model does not appear to 
be able to account for the metastability of the supercooled liquid to any significant degree, 
so providing tacit support for the idea that some sort of competition between multiple stable 
configurations is required to properly describe the supercooled state. We shall return to this 
point in the conclusion. In this paper we shall explore the consequences of the minimal model 
on the liquid and solid states. Our primary interest is to understand how the geometry of 
the favoured local structure (FLS) determines the nature of the freezing transition. 

To address the relationship between the liquid structure and freezing, we shall use a 2D 
lattice model, the Favoured Local Structure (FLS) model, which allows us to directly select 
the geometry of a favoured local structure from among the possible geometries^. The benefit 
of the FLS model is that we do not need to wrestle with the non-trivial task of trying to 
engineer a desired local geometry through particle interactions. This task has been the 
subject of a number of recent studied^. A specific challenge is that, for small clusters 
of particles interacting by non-directional short range potentials, the energy proves to be 
rather insensitive to cluster geometry. This has been addressed by either adding multiple 
components (up to 8 in an 8 particle clusteiP) or by exploiting long-range interaction^ 
to select a desired cluster geometry. By assigning energy explicitly based on the cluster 
geometry itself, as we do in this paper, we effectively set aside the problem of engineering a 
set of particle interactions to achieve this stability so that we can focus on the consequences 
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- thermodynamic and kinetic - of the symmetry of a particular favoured local structure. 

Previously, we have identified the crystal ground states for each possible choice of FLS 3 
and demonstrated that the decrease in liquid entropy with decreasing energy is smaller for 
liquids based on a low symmetry FLS as compared with those based on high symmetry local 
structures 6 . In this paper, we present an account of the freezing transitions described by the 
FLS model. In doing so we shall demonstrate how the microcanonical treatment of the phase 
transitions allows us to clearly see the dual role of the local structure in both establishing 
the energy of the crystal state and the entropy cost of local ordering in the liquid state. 

II. THE FLS MODEL 

We consider a 2D triangular lattice of Ising spins (each site has a single degree of freedom 
named spin value, which can take two values called up or down). Readers should note that 
this model is identical to one in which two different species A and B are distributed on the 
lattice. The local environment of a given site is defined as the spin state of its 6 nearest 
neighbours. In the set of 2 6 possible environments, there are 8 local structures that cannot 
be inter-converted by rotations or spin inversion. These structures are sketched in the insets 
on Figure [2} The labeling convention we adopted to design each local structure goes as 
follow: the first digit is the number of down spin (dark sites on the pictures), and the second 
one is the length of the longest sequence of up spins in the structure (when there is only or 
1 down site, this second digit is unnecessary, and associated structures will be named simply 
{0} and {1} respectively). The {32} structure has no plane symmetry and is therefore chiral. 
Here we will consider only the case when the enantiomers are not distinguished. The chiral 
case {32c}, where we favour only one enantiomer, involves a number of interesting subtleties 
that will be addressed in a separate papei^. 

Among these 8 distinct local structures, we select one that will be designated as the 
favoured local structure (FLS). Each site whose local environment is the FLS (to a rotation) 
is assigned an energy of — 1 ; all the others have an energy of 0. Note that this energy is 
independent of the spin on the given site, depending only on the spin arrangement of the 6 
neighbouring sites. 

We evaluate the equilibrium properties for each FSL model using Monte Carlo sampling. 
We consider a system in the canonical ensemble, at temperature T (in this work we will 
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take fee = 1). For high temperatures we have used the standard Metropolis algorithm and, 
at low temperatures, the equivalent rejection- free method due to Bortz et al 8 in which sites 
are organized in terms of the energy change associated with the spin flip. 

We have previously 3 identified the ground states (all crystalline) for each choice of FLS. 
To assist the discussion of the various freezing transitions, we present the crystal phases 
for the different FLS's in Figure [T] along with the energy per site E Q and the number g of 
distinct realizations of a FLS on the lattice. Note that three FLS's have two distinct crystals 
with the same energy. 

III. LIQUID-CRYSTAL PHASE TRANSITIONS AND THE MICROCANONICAL 
MAXWELL CONSTRUCTION 

The temperature dependence of the average energies on heating and cooling are plotted 
in Figure [2] for each FLS. We note that the choice of FLS can result in either first order 
transitions ({1}, {24} and {33}), continuous transitions ({0}, {22} and {31}), or, in the 
case of {32} and {23}, multiple transitions. In the cases where two degenerate polymorphs 
exist we find that only one of the two is observed to form in the case of {23} and {24}. 
The observed crystal is labeled a in Figure [TJ In the case of {31}, the polymorphs are 
spin-inversion symmetric and either can form spontaneously from the liquid. 

In addition to the variety of freezing transitions, the equilibrium liquid states exhibit a 
considerable variation in their heat capacities Cy = and the maximum concentration of 
FLS that the liquid can accumulate before freezing. In a previous paper 6 we demonstrated 
that the influence of the choice of FLS on the liquid structure could be attributed to the 
difference in the entropy cost per FLS. Contributions to this entropy cost per FLS were 
found to come from both the symmetry of the favoured local structure itself (the higher the 
symmetry, the greater the entropy cost) and from the entropy of aggregates of FLS's. 

Not only does the decrease in symmetry of the FLS allow the liquid to accumulate more 
local order before freezing, but it also tends to result in large unit cells in the crystaP. This 
coincidence of effects means that large unit cells in the crystal phase will often be associated 
with a liquid that can accumulate a significant amount of local structure (i.e. 'pre-order') 
before freezing and help explain why we see little difference in the rate of crystallization of 
the different systems, even though the unit cell varies across a factor of 20. To appreciate 
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the extent of this pre-ordering, consider the case of the {32} FLS which freezes into our 
most complicated crystal structure (see Figure [I] This freezing involves multiple transitions 
(as we shall describe below) but if we select T = 0.6, a temperature that lies just above the 
first of the ordering transitions, we find that the local ordering has already reached roughly 
75% of the final crystal order. In Figure [3] we show an example configuration in the {32} 
liquid at T = 0.6 in which the entangled fragments of the crystal order are clearly evident. 

Given the established connection between the structure of the FLS and the dependence 
of entropy S(E) on E in the liquid, it seems natural to consider the freezing transitions 
from the perspective of the microcanonical ensemble. In this approach, coexistence arises 
from maximising the entropy at a fixed value of the energy. We can always recover the 
temperature from the relation 

dE T K ; 

Since both the energy and the entropy are bounded, this approach allows us to embrace 
the entire temperature range in one plot. In the thermodynamic limit, this construction is 
completely equivalent to canonical-ensemble constructions. 

We can calculate the entropy by integrating 1/T over the energy, i.e. 

pE jp/ 

S(E) = s{e ' ) + ] e ,tW) (2) 

where T(E) is the temperature associated to the average energy E in Figure [5J To carry out 
this calculation we need to know the value of the entropy S(Ei) at the value of the energy 
Ei from which the integration is carried out. The liquid state is associated to the high 
temperature limit. At T — » oc, the entropy per spin is = ln(2) (the factor 2 reflecting 
the 2 possible states of any individual site) and the energy is the average concentration of 
FLS's in a random state of the system, i.e. 

Eoo = -| (3) 

where g is the degeneracy of the FLS and is provided in Figure [T] For the solid entropy, 
the ground state energy E (also provided in Figure [j]) provides the reference energy with 
an entropy S(E ) = 0. By using Eq. [2] and the respective references states for the liquid 
and solid, we have calculated the liquid and solid entropies versus energy for each choice of 
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FLS. Any state with an entropy below the convex hull of the S(E) curve corresponds to a 
metastable state. The common tangent to liquid and crystal entropies therefore corresponds 
to the thermodynamic freezing transition, and its slope is 1/7/. We have constructed this 
common tangent with the help of fitted curves for the entropy (power law at low temperature, 
polynomial for the liquid part). This construction is presented on Figure [4] for each choice 
of FLS. 

The transition temperature for the continuous transitions are taken from the peak in the 
heat capacity. The values of Tf for each of the choices of FLS are presented in Table [I] 
along the values of Tf as determined by the Maxwell construction. The common tangent 
construction is strictly inapplicable to second order transitions. In practical terms, any 
attempt to identify distinct solid and liquid entropy curves for continuous transitions runs 
into the problem of sensitivity to the specific choice of energy intervals used to identify the 
two phases. That said, the Maxwell construction, applied to the continuous transitions, does 
provide a quite reasonable estimate of the transition temperature. 





Continuous 


First order 


Multi-step 


FLS 
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{31} 


{22} 


{33} 


{24} 


{1} 


{32} 


{23} 


Tf 


1.43 


0.85 


0.96 


0.69 


0.57 


0.53 


0.34 


0.53 
















0.50 


0.57 
















0.56 


0.64 


Tf, maxwell 


1.43 


0.87 


0.97 


0.70 


0.57 


0.54 


0.52 


0.60 



TABLE I: The freezing points of the various FLS systems. In the case of a continuous transition 
Tf is determined as the position of the maximum of the heat capacity. In the case of a first 
order transition, Tf is identified as the midpoint between the transition temperatures on hear and 
cooling. In case of multiple transitions, all the individual transition temperatures are presented. 
^/, maxwell was obtained via the common tangent construction as shown in Figure [4j 

The microcanonical representation of the freezing transition allows for the most transpar- 
ent connection to be made between the immediate consequences of the choice of FLS and 
the transition temperature. To illustrate this connection, we shall consider the following 
simple treatment of the freezing transition. Let the liquid entropy by given by, 
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S(E) = - A(E 00 - Ef (4) 

As explained in the Appendix, Eq. [| represents the first term in an expansion of the liquid 
in terms of (E^ — E) and the constant A can be evaluated exactly for this expansion. These 
calculated values of A are provided in ref. 6 . The second assumption is to neglect the entropy 
of the crystal, and assume that its entropy at coexistence is zero. This assumption is justified 
by the observation on Figure [4] that the entropy range for the crystal is indeed very small 
as compared to the liquid's. A more accurate treatment of the crystal entropy is provided 
by a single excitation model discussed in the Appendix. With these two assumptions, the 
common tangent construction is expressed by the relation 

2A( Eoo — E) — S ~-f^ E ~ E? (5) 
Solving for E, the value of the liquid energy at coexistence, we have 

E = E + VA 2 - S~/A (6) 

where A = E^ — E Q . The theoretical estimate of the transition temperature 7^ /ieor2/ 5 given 
by the inverse slope of the common tangent, is 

rjitheory 1 /y\ 

f ~ 2A(A - ^A 2 - SUA) 
As is evident from Eqs. [6]and[7j the freezing transition is determined (within this approxi- 
mate treatment) by just two parameters: the crystal lattice energy, E^ — E Q) and the liquid 
entropy factor A. The latter factor contains the entropy cost of structure in the liquid : 
a large value of A ) generally associated with a high symmetry FLS, corresponds to a local 
structure that incurs a large entropy cost. In Fig. [5] we plot the values of T^ heory calculated 
using Eq. [7] with the values of the freezing temperature from Table [1} Given the simplicity of 
the treatment, the theoretical estimates of the freezing temperatures are quite reasonable. 
The notable exception to this success is the FLS {1} where the theory considerably over- 
estimates transition temperature. The {1} liquid actually exhibits a substantial decrease 
in energy prior to freezing, which is not correctly taken into account by the factor A ) as 
previously noticecP, which explains this deviation. While this approximate treatment is too 
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simple to establish whether a transition is first- or second-order, we note that it does ap- 
pear to get the magnitude of the transition temperature roughly right, even for continuous 
transitions. 

IV. FREEZING IN MULTIPLE STEPS 

The presence of multiple transitions for the {23} and {32} FLS's is clearly seen when we 
examine the E(T) curves with an expanded temperature scale about the transition region, 
as shown in Figure [6| In both cases, the freezing transition is made up of three distinct 
transitions. This is unusual for crystallization which is typically viewed as the archetype of 
a highly cooperative transition and hence unlikely to decouple in this way. In this Section we 
shall try to identify the intermediate phases associated with the step- wise freezing processes. 

In the case of the {23} FLS, the multiple transition scenario appears to arise from the 
existence of the two degenerate polymorphs (see Figure [T]). If we identify the three transition 
temperature as 7\, T 2 and T 3 , as indicated on Figure [6], it appears that the transition on 
cooling through T 2 corresponds to the formation of a polycrystalline mixture of the two 
polymorphs. Further cooling sees the formation at T 3 of the single phase of the preferred 
crystal phase. The nature of the continuous transition at T\ remains a puzzle. 

Turning to the transitions in the {32} system , the picture is quite different. Instead of 
different crystal forms, the sequence of transitions correspond to a step-by-step symmetry 
breaking. Transition 1 marks the appearance of an orientation, at which the symmetry of the 
liquid is reduced from "isotropic" to a liquid with 3-fold rotational symmetry. On cooling, 
transition 2 is characterised by the loss of this rotational symmetry as unique orientation is 
selected. Finally, transition 3 breaks the reflection symmetry and the translational invariance 
in one direction, to give the low T crystalline state. Transitions 2 and 3 seem to be first 
order, though the small energy gap at transition 3 permits oscillation between phases in 
our finite size systems. This sequence of symmetry breakings is reminiscent of the liquid- 
nematic-smectic sequence of transitions observed in some molecular liquids. 
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V. CONCLUSION 



In this paper we have examined the variety of freezing transitions exhibited by a model 
liquid characterised by a favoured local structure. We found examples of ordering transitions 
that were first order, continuous and multi-stepped. All liquids were found to order with 
little or, in the case of the continuous transitions, no evidence of a tendency to supercool, in 
spite of the involvement of some complex crystal phases. The intuitive notion that large unit 
cell crystals are more 'difficult' to crystallize is challenged here with the demonstration that 
the very same features that produce such crystals also ensure that a significant amount of 
the low energy structures will be present in the liquid in advance of the freezing transition. 

The microcanonical representation of the freezing transition was presented as the most 
transparent way of expressing the influence of the local structure on the crystal (through 
the groundstate energy E Q ) and the liquid (via the factor A related to the entropy cost per 
FLS) on the freezing transition. A simple theoretical treatment that retained only E^ — E Q 
and A was found to provide a reasonable estimate of the transition temperature for most of 
the liquids considered. 

The FLS model involves considerable simplifications. Some, like the restriction to just 
one FLS, can be removed easily. Others, such as the neglect of the fluctuation of the local 
structures away from their ideal form cannot easily be included due to the restriction of 
the lattice. This work poses the question - is there a useful way of expressing, at least 
approximately, an effective Hamiltonian for interacting particles with continuous degrees of 
freedom in terms of explicit favoured local structures. 

By construction, the crystalline phase of our model is the densest packing of FLS. For all 
FLS's other than {0}, these packings don't completely tile the space, a consequence of the 
fact that these FLS's are frustrated. A large body of work on the notion of geometrical frus- 
tration, introduced by Frank in 1952 commentP, rests on the proposition that the presence 
of a frustrated FLS can be used to account for the metastability of the supercooled liquid. 
Based on our results for the FLS model, frustration of an FLS alone is generally insufficient 
to stabilize the supercooled liquid. The competition between two FLS's: one that can pro- 
duce a low energy crystal while the other FLS, although more stable, can only order in a 
higher energy structure, looks like the simplest means of generating a glass forming liquid. 

While acknowledging the limitations of the FLS model, we believe that it does provide an 
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excellent tool for building up our intuition on how local structure influences the properties 
of condensed matter. A report on the extension to 3D is currently in preparation^. The 
role of competing local structures in stabilizing the supercooled liquid and influencing the 
kinetics of crystallization would involve a straight forward extension of the model presented 
here. By building the model around the idea of local stable structures, rather than having 
to discover them as a consequence of particle interactions, we believe that the FLS approach 
offers a clearer language with which to ask questions about structure in condensed phases - 
both the influence of local structure and the nature of global ordering. 
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APPENDIX 

In this Appendix we describe an alternative to the numerical fits used for the common 
tangent construction, which gives analytic values of the Taylor expansion coefficients around 
asymptotic high- and low-temperature regimes. We formerly derive the first order term of 
the expansioiM Here we generalize these results to any order using correlation functions 
that can be computed by exact enumeration. We also present a route to approximate the 
thermodynamic quantities at low temperature, though only to leading order. 
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A. The Energy of the Low Temperature Crystal 

We shall treat the energy of the crystal phase using a low temperature approximation 
in which we only include the (uncorrelated) single spin excitations from the groundstate. 
First, we identify the energy spectrum {n i) Ej) of the ground state, as introduced in Table 
in the example of system {22}. It classifies the single-spin defects that can occur in the 
crystal according to the energy of the corresponding state Si and the proportion of site that 
can host such a defect r^. 



State i 


• •••• 

)«0»OI 

• •••• 

>o«o»c 
••••• 


• •••• 

)#0«OI 

••••• 

>o»o»c 
••••• 


• •••• 

)#0#OI 

• •••• 
> ••• 

• •••• 


Energy cost 


(ground state) 


4 


6 


Concentration n$ 




3/4 


1/1 



TABLE II: Sketch of the energy spectrum of system {22} at low temperature. For this system, two 
kinds of elementary (single spin) defects can occur in the crystal, examples of which are presented 
schematically (the ground state is also reminded). The energy spectrum of a crystal is determined 
by the energy costs e% of each kind of defect, and the respective concentration of sites that can host 
such a defect n^. 

In the canonical ensemble at temperature T, each site of type % has a probability: 

to be in its excited state, if we negligate the interaction between such defects. This gives 
the following approximate expression for the energy : 

E ( T ) = E ° + T,^ 1 + e-l/r (8) 

i 

All other thermodynamic quantities derive from E(T). However the computation of S(T), 
for example, involves complicated analytic forms, and it is much easier to compute it by 
numerical integration of Eqj8] using Eq. [2j to arbitrary precision. 

A quantitative comparison between the results of simulations and this model is presented 
on Figure[7J for system {22}. We find that the accordance is quantitative up to a temperature 
around T = 0.6. For the other systems, no qualitative difference was observed, though 
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the quantitative values in the energy spectrum vary. This simple model gives a satisfying 
description of the crystal state, the physics of which is therefore quite simple. This method 
gives easily the leading order of the low-temperature expansion. Finding the corrections to 
this expansion would involve correlation between defects. No difficulty is expected to arise 
from that, though the computations would be rather tedious. 



B. Exact high-temperature expansion of the FLS model at all orders 

In this section, we present an adapted version of the classical high-temperature expansions 
of spin models adapted to the FLS model. It gives the exact values of the coefficients of the 
Taylor expansion of the thermodynamic functions around f3 = 1/T = 0. These coefficients 
are connected correlation functions at infinite temperature, which may be computed exactly. 
This reasoning is an extension of our work in 6 , where we presented the route to compute 
the first non-trivial coefficient. 



1. Definitions : 

To make the link with other spin models, we first introduce the Hamiltonian (energy of 
a configuration) : 

= (9) 
where A is the lattice, C the spin configuration, and the local energy variables are : 

{ — 1 if i is in the FLS in configuration C 
(10) 
otherwise 

Each site of the lattice has z neighbours [z is the coordinence of the lattice). We assume 
invariance under translations, which impose to use periodic boundary conditions if the sys- 
tem is of finite size. For the triangular lattice studied here z = 6. Note that the FLS model 
could be adapted to other regular lattices ; a review of this model on the 3D FCC lattice is 
in preparation. 

The partition function of the system is defined as : 

Z(f3) = J2eM-PH(C)) (11) 
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where f3 — 1/T is the inverse temperature. All thermodynamic quantities derive from Z 
and its derivatives ; for example : 

F(f3)= -JlogZGS) (12) 

m= -zb) 9 -T (13) 

S(P)= p(E-F) (14) 

are respectively the free energy, the average energy and the entropy of the system. 

We define an average (-) as the uniform average over all configurations, i.e. at /3 = 0. 
Note that for any site z, we have : 

(e i ) = E oo = -^ (15) 

where g is the degeneracy of the FLS and is the average energy per site at infinite 
temperature. 



2. Derivation 



Now that the notations are set, we are ready to derive an exact Taylor expansion, at any 
order, of the thermodynamic functions in the FLS model around /3 — 0. 

Working on the partition function Z(/3), a few usual tricks allow us to rewrite: 

Z(J3) =£ c exp(-/3#(C)) (16) 

= E C ELeA ex P(/^) ( 17 ) 
= Ec a e A(l + vet) (18) 

with v — e 13 — 1. Note that v w /3 at high temperature. These simplifications use the fact 
that q takes only two values, and 1. 

Now we can expand this product exactly : 

na + ^) = $> |A| ]> (19) 

ieA ACA ieA 

The meaning of this expansion is the following : A is the subset of right-side factors se- 
lected in the expansion of the product, and \A\ is the cardinal of A. If we now sum over 
configurations C : 

Z{P)= EcUieA^ + ^i) (20) 
= 2^EACA^<rLeWo (21) 

13 



where TV is the total number of sites. We can rearrange this sum according to \A\ : 

/cGNACA, \A\=k ieA 

So finally we get the high temperature explicit, exact expansion : 



(22) 



(23) 



A convenient rewriting of this equation in terms of the free energy leads us to introduce 
the connected correlation functions: 



(3F((3) = - log Z((3) = -iVlog(2) - E •••**> 



k>0 ix^—^ik 



(24) 



where : 



fe>o = (e<>o 
( 6 i 6 j)o 7 = ( 6 * e j)o - (q)o(^-)o 
( e i e j e k)o = (ti£j£k)o — (ci)o(ejek)o ~ ( e j)o( e i e k)o — {ck)o(€i€j)o + 2(e^)o(ej)o(e^)o 

For the explicit computation, free energy representation is more convenient since con- 
nected correlation functions vanish as soon as a site is "far" (not overlapping) from the 
others. Thus every term in the expansion is extensive in the size of the system. In order 
to compute explicitly the connected correlation function at order fc, we need to enumerate 
all possible ways to put an FLS in each of the k sites ii . . . i^. It will be zero if they can 
be separated in two non-overlapping set of structures, thus we need only to enumerate a 
finite number of positions i\ . . . i^. The quantity (e^e^ . . . e ik )o can be interpreted simply 
as the probability that all sites i\ . . .i^ lie simultaneously in the FLS, in a random spin 
configuration. 

3. Computation of the first terms 



The first terms read : 

PFifi) = -AHog(2) - + S E< e * e ^ + 1/3 E te £ i e *>o + 0(v A ) 



(25) 



14 



The truncature to the order v 2 would give the Gaussian approximation to the density of 
states, as iiP. Now the interest of this new, more formal approach, is to be able to increase 
the precision and involve higher-order correlation functions. 

All that is left is to compute explicitly these connected correlation functions, which need 
only to test a finite (but possibly big) number of configurations. This is the only model- 
specific part of this derivation. Using enumerating routines, we could quite easily get the is 3 
term for FCC and triangular lattices. Finally, to obtain an expansion of the free energy in 
powers of /3 instead of z/, we can substitute the expansion of v in terms of /3, i.e. 




(26) 



fc>0 



into Eq. [25] 
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FIG. 1: The groundstates for the 8 distinct FLS as defined in the text. The geometrical multiplicity 
g of the FLS, number of sites in the unit cell Z and the energy per site Eq in the groundstate are 
indicated. For the {23}, {24} and {31} FLS's we find two degenerate polymorphs. 
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FIG. 2: The average energy per site as a function of temperature for each FLS. The measurements 
were made in an hysteresis cycle of temperature -cooling (in blue) then heating (in red). These 
curves were obtained with a 60 x 60 lattice and a 'cooling rate" of 2.10 6 steps of the rejection- free 
algorithm per site and per unit of temperature. 
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FIG. 3: Typical equilibrium spin configuration of the system for FLS {32}, at temperature 
T — 0.6, just above the first transition. This liquid state is highly structured. 
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FIG. 4: The entropy S(E) for the crystal (red) and liquid (blue) for each choice of the FLS. The 
common tangent is indicated as a dashed line in each case. 
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FIG. 5: A scatter plot of the values of Tf from simulation against the value of T^ heory from Eq. 
The straight line with slope one indicates the case where Tf = 7^ heory m 
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FIG. 7: The low temperature model described in equation [8j compared to simulated results for 
system {22}. The agreement is excellent at low temperature, as can be seen on the right figure, 
which is E(T) — Eq plotted in logarithmic scale. This model correctly identifies leading-order 
behaviour at low temperature. 
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